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Abstract 

An effective field theory approach is developed for calculating the thermodynamic 
properties of a field theory at high temperature T and weak coupling g. The effec- 
tive theory is the 3-dimensional field theory obtained by dimensional reduction to the 
bosonic zero- frequency modes. The parameters of the effective theory can be calculated 
as perturbation series in the running coupling constant g^iT). The free energy is sep- 
arated into the contributions from the momentum scales T and gT, respectively. The 
first term can be written as a perturbation series in g'^iT). If all forces are screened 
at the scale gT, the second term can be calculated as a perturbation series in g{T) 
beginning at order g^ . The parameters of the effective theory satisfy renormalization 
group equations that can be used to sum up leading logarithms of T/{gT). We apply 
this method to a massless scalar field with a <I>'^ interaction, calculating the free energy 
to order g^ log (7 and the screening mass to order g^ logg. 



I. Introduction 



When matter is subjected to sufficiently extreme temperature or density, both quantum 
and relativistic effects become important. Such conditions arise in astrophysics and in cos- 
mology and it may be possible to create them experimentally using heavy-ion collisions. A 
system at such an extreme temperature or density is most appropriately described by quan- 
tum field theory. The thermodynamic functions that describe the bulk equilibrium properties 
of such a system are given by the free energy density and its derivatives. If the temperature 
T is high enough that all masses can be neglected, the free energy density depends only on 
T and on the coupling constants of the field theory. 

In recent years, there have been significant advances in the perturbative calculation of 
the free energy for high temperature field theories. The free energy for a massless scalar field 
with a interaction was computed to order g'^ by Frenkel, Saa, and Taylor in 1992, and 
the order-(7^ correction was recently calculated by Parwani and Singh 0. The free energy for 
high temperature QED was calculated to order by Coriano and Parwani , and extended 
to order by Parwani . The free energy for a quark-gluon plasma in the high temperature 
limit was recently calculated to order by Arnold and Zhai [Q. The importance of these 
calculations goes far beyond simply determining one more term in the perturbation series. 
The leading term in the perturbation series is just the free energy of an ideal gas. The order- 
g"^ correction takes into account interactions between particles in the ideal gas. At order 
g^, there is a qualitatively new contribution to the free energy. If the temperature is large 
compared to the masses of the particles, the force mediated by the exchange of a particle is 
long-range compared to the typical separation of the particles, which is of order 1/T. The 
system therefore behaves like a plasma, screening the long-range interaction beyond the scale 
l/{gT). It is this screening that is responsible for the g^ term in the free energy density. 
Because of renormalization effects, the corrections of order g^ and g^ are also important. 
The coupling constant g{fi) depends on an arbitrary renormalization scale fi. The resulting 
ambiguity in the leading nontrivial term in a perturbative expansion can only be reduced 
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by a next-to-leading order calculation. While it may be clear on physical grounds that the 
scale IX should be of order T, the difference between the choices g{T) and g{27iT) can be of 
great practical significance. Thus a calculation to order is needed in order to determine 
the appropriate scale /i in the ordeT-g^ correction to the ideal gas term, while a calculation 
to order g^ is required in order to determine the scale in the order-(yf^ plasma term. 

In the case of QCD, there is also a qualitatively new effect that arises at order g^. As 
pointed out by Linde 0, in 1979, the loop expansion for the free energy breaks down 
at this order in g. Up to order g^, the free energy can be calculated using a resummation 
of perturbation theory that takes into account the screening of the chromoelectric force at 
distances of order l/{gT). However the chromomagnetic force is not screened at the scale 
gT, and this causes a breakdown in the resummed perturbation expansion at order g^. This 
breakdown has been widely interpreted as implying that nonperturbative effects enter in at 
this order, rendering the perturbation series meaningless beyond order g^. This longstanding 
problem was recently solved by constructing a sequence of two effective field theories that 
are equivalent to thermal QCD over successively longer length scales |^. The first effective 
theory reproduces the static gauge-invariant correlators of thermal QCD at distances of order 
1 / {gT) or larger, while the second effective theory reproduces the correlators at distances 
of order l/{g'^T) or larger. Using this construction, the free energy can be separated into 
contributions from the momentum scales T, gT, and g^T, with well-defined weak coupling 
expansions that begin at order g^, g^, and g^, respectively. The contributions from the 
scales T and gT can be computed using perturbative methods. For the contribution from 
the scale g'^T, the coefficients in the weak-coupling expansion can be calculated using lattice 
simulations of pure-gauge QCD in 3 Euclidean dimensions. This result demonstrates the 
power of the effective-field-theory approach. This power can also be brought to bear on 
other fundamental problems in thermal field theory. For example, it has also been used to 
determine the correct asymptotic behavior of the correlator of Polyakov loop operators in 
high temperature QCD 
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Effective field tfieory provides an effective metliod for unravelling the effects of the various 
momentum scales that arise in field theory at high temperature. In addition to providing 
insight into the qualitative behavior of the theory, effective field theory can also be used 
to streamline perturbative calculations, such as those in Refs. In this paper, we 

use effective field theory to develop a practical method for calculating the thermodynamic 
functions of a field theory at high temperature T and weak coupling g. To illustrate the 
method, we apply it to a massless scalar field theory with a interaction. In section II, we 
explain how the effective field theory that describes distance scales of order 1 / {gT) or larger 
is related to dimensional reduction to the bosonic zero-frequency modes. In Section III, we 
show how the parameters of the effective theory can be obtained by matching perturbative 
calculations in the effective theory and in the full theory. In Section IV, we use the effective 
theory for the massless $^ theory to calculate the free energy to order g^ and the screening 
mass to order g'^. The accuracy of these calculations is improved in Section V to order g^ log g 
for the free energy and to order g^ log g for the screening mass by using renormalization 
group equations for the parameters of the effective theory. We conclude in Section VI with a 
discussion of the application of our effective-field-theory approach to gauge theories and with 
a brief comparison with related work. In appendices A and B, we collect all the necessary 
formulas for sum-integrals in the full theory and for integrals in the 3-dimensional effective 
theory. In appendix C, we derive the renormalization group equations for the parameters of 
the effective theory. 



II. Dimensional Reduction 

In the limit of high temperature T, the static correlation functions of a field theory in 
(3+1) dimensions can be reproduced at long distances i? ^ 1/T by an effective field theory 
in 3 dimensions. This idea, which is called "dimensional reduction", has a long history 
1^, 10, 11, 12]. It has provided insight into the qualitative behavior of the field theory at 



high temperature, but it has never been fully exploited for quantitative calculations. 



Dimensional reduction is based on the fact that static correlation functions for a field 
theory in thermal equilibrium can be expressed in terms of Euclidean functional integrals. 
The partition function Z is defined by 

Z{T) = trace (e"^^) , (1) 

where H is the hamiltonian operator and P = 1/T. The operator e~^^ is the evolution 
operator that evolves a state from time t = to the imaginary time t = —i(3. It can be 
represented as a functional integral over fields $(x, t) defined on the time interval from to 
—i(3. It is natural to change variables from t to the imaginary time r = it. The partition 
function ([l|) is then given by the Euclidean functional integral 

Z{T) = J V^{^,t) exp J d^x , (2) 

where C is the negative of the lagrangian density for the (3+l)-dimensional theory with the 
time t analytically continued to —ir. The trace in is implemented by imposing boundary 
conditions on the fields: 

$(x,r = /5) = ±<l>(x,r = 0), (3) 

where the plus sign holds for bosonic fields and the minus sign for fermionic fields. The 
correlator of two operators O{0) and 0(R) is obtained by averaging their product over 
fields $(x, r) with the exponential weighting factor in (|^). 

Because of the periodicity conditions (|^) on the fields, they can be decomposed into 
Fourier modes in r, with Matsubara frequencies a;„ = 2mTT for bosons and uJn = (2n + l)7rT 
for fermions. The contribution to a correlator from the exchange of a Fourier mode with 
frequency LOn falls off at large R like exp{~\uJn\R) ■ Thus the only modes whose contributions 
do not fall off exponentially at distances greater than 1/T are the n = modes of the 
bosons. This suggests the strategy of integrating out the fermionic modes and the nonzero 
modes of the bosons to get an effective theory for the bosonic zero modes. This process 
is called "dimensional reduction" . It results in a 3-dimensional Euclidean field theory with 
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bosonic fields only which reproduces the static correlators of the original theory at distances 
R > l/T. 

Constructing the dimensionally-reduced effective theory by actually integrating out de- 
grees of freedom is cumbersome beyond leading order in the coupling constant. Once the 
appropriate 3-dimensional fields and their symmetries have been identified, a better strategy 



is to use the methods of "effective field theory" |13|. One writes down the most general 
lagrangian £eff for the 3-dimensional fields that respects the symmetries. This effective la- 
grangian has infinitely many parameters, but they are not arbitrary. By computing static 
correlators in the full theory, computing the corresponding correlators in the effective theory, 
and demanding that they agree at distances R ^ 1/7", one can determine the parameters 
of Ccs in terms of T and the parameters of the original theory. Note that this matching 
procedure does not necessarily require the explicit determination of the relation between the 
fields in the effective theory and the fundamental fields. 

The construction of the 3-dimensional effective theory is complicated by ultraviolet di- 
vergences. The ultraviolet divergences associated with the original 4-dimensional theory are 
removed by the standard renormalization procedure. However the 3-dimensional effective 
theory also has ultraviolet divergences whose origin can be traced to integrating out the 
nonzero-frequency modes in the full theory. They must be regularized by introducing an 
ultraviolet cutoff A. The parameters in £eff must therefore depend on A in such a way as 
to cancel the A-dependence of the regularized loop integrals in the effective theory. The 
natural scale for the ultraviolet cutoff A is of order T, since this is the scale at which the 
corresponding integrals are cut off by the nonzero modes in the full theory. It is useful, 
however, to keep the cutoff A independent of T, so that the only dependence on T in the 
effective theory comes from the parameters in £efr- From the point of view of the full theory, 
the ultraviolet cutoff A of the effective theory plays the role of an arbitrary factorization 
scale that is introduced to separate the momentum scale T from lower momentum scales, 
such as gT, which can be described within the effective theory. 
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The ultraviolet divergences of the effective theory include power ultraviolet divergences 
of the form A^, p = 1,2,3,..., and logarithmic ultraviolet divergences of the form log(A/m), 
where m is a mass scale in the effective theory. The power divergences and the logarithmic 
divergences are quite different in character. The coefficients of the power divergences depend 
on the regularization procedure and are therefore simply regularization artifacts. In particu- 
lar, the coefficient of in some observable calculated in the effective theory is independent 
of the coefficient of for that observable in the full theory. Therefore power divergences 
of the form from loop integrals must be completely cancelled by terms proportional to 
A^ in the parameters of C^s- In contrast to the power divergences, logarithmic ultraviolet 
divergences have coefficients that are independent of the regularization procedure. The rea- 
son for this is that a logarithmic ultraviolet divergence of the form log(A/m) from a loop 
integral must match onto a log(T/A) term in one of the parameters of C^s in order for the 
A-dependence to cancel. Thus the logarithmic divergences in the effective theory are related 
to logarithms of T in the full theory and therefore have real physical significance. Because of 
the unphysical character of power ultraviolet divergences, it is convenient to use a regulariza- 
tion procedure for the effective theory in which power ultraviolet divergences are subtracted 
and the remaining logarithmically divergent integrals are cut off at the scale A. With such 
a regularization procedure, only the logarithmic divergences need to be cancelled by the pa- 
rameters of the effective theory. The subtraction of power divergences can be implemented 
with any cutoff procedure, including a momentum cutoff or lattice regularization. For pertur- 
bative calculations, a particularly convenient cutoff procedure is dimensional regularization, 
in which momentum integrals are analytically continued to 3 — 2e spacial dimensions. Power 
divergences are automatically subtracted with this method, because integrals without any 
momentum scale are by definition in dimensional regularization. The remaining logarith- 
mic ultraviolet divergences appear as poles in e, and the renormalization procedure can be 
completed by subtracting these poles. If A is the momentum scale introduced by dimen- 
sional regularization, then this "minimal subtraction" procedure is equivalent to cutting off 
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the logarithmic ultraviolet divergences at a momentum scale of order A. 

We turn now to the specific case of a massless scalar field with a interaction. In the 
partition function (0), the Euclidean lagrangian density is 

L = \ [dM + \ (V$)^ + (4) 

The 3-dimensional effective field theory obtained by dimensional reduction describes a scalar 
field 0(x) that can be approximately identified with the zero-frequency mode of the field in 
the original theory: 

Vt / dr $(x,r) ^ 0(x). (5) 

JO 

The symmetries of the effective theory are rotational symmetry and the discrete symmetry 
0(x) — 0(x), which follows from the symmetry $ — $ of the fundamental lagrangian 
(^). The effective lagrangian has the general form 

£eff = ^(V0)^ + im2(A)0^ + ^A(A)0^ + 5£, (6) 

where 5C includes all other local terms that are consistent with the symmetries. The param- 
eters m^(A) and A(A) and the infinitely many parameters in depend on the ultraviolet 
cutoff A, the temperature T, and the coupling constant g"^. The partition function can 
be expressed as a Euclidean functional integral over the 3-dimensional field 0(x): 

Z{T) = e-^(^) ^ 1^^^ l)0(x) exp J d'x C,s) ■ (7) 

The parameter /(A) in the exponential prefactor is the coefficient of the unit operator in the 
effective lagrangian, which was omitted from £efr in (H)- In addition to depending on g"^ and 
T, it also depends on A in such a way as to cancel the A-dependence of the functional integral 
in (^. The correlator of operators O{0) and 0{TV) in the full theory can be calculated at 
long distances R ^ 1/T by identifying the corresponding operators in the effective theory 
and averaging their product over fields 0(x) with the exponential weighting factor in (|^. 

Renormalization theory implies that correlators at long distances R ^ 1/T can be re- 
produced to any desired accuracy by adding sufficiently many operators to the effective 
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lagrangian and tuning their coefficients with sufficient accuracy as functions of g'^, T, and A 
|r3| . With the three terms in the lagrangian that are given exphcitly in long-distance 
correlators can only be reproduced with limited accuracy. If we include the operator 0^, the 
field theory becomes renormalizable, but this renormalizable theory is still only accurate up 
to a finite order in the coupling constant g. This has been interpreted as a breakdown of 



dimensional reduction |12|, but the correct interpretation is simply that nonrenormalizable 
operators must also be included in Ces in order to extend the accuracy to higher order in 
g. In general, one must include all operators that are invariant under rotations and under 
the symmetry —0. The resulting field theory is nonrenormalizable and has infinitely 
many parameters. These parameters, however, are not arbitrary, but are determined as func- 
tions of g^, T, and A by the condition that the effective theory reproduce the long distance 
behavior of the original theory. 

It is easy to determine the magnitude of the coefficient of a general operator in the 
effective lagrangian. From the kinetic term (V0)^ in (^), we see that the field should 
be assigned a scaling dimension of 1/2. The operators given explicitly in then have 
dimensions 3, 1, and 2. If we use a renormalization procedure in which power ultraviolet 
divergences are subtracted, then by dimensional analysis, an operator of dimension d in the 
effective lagrangian must have a coefficient that is proportional to T^^'^. It remains only to 
determine its order in g. The operator 0^ is generated at tree level from the g^^"^ term in the 
original lagrangian, and therefore has a coefficient proportional to g^. All other interaction 
terms in Ces arise from loop diagrams in the effective theory. Operators with 2n powers of 
are generated by 1-loop diagrams with n 4-point interactions and therefore have coefficients of 
order g'^"'. Thus an operator with the schematic structure V^™'0^"' has dimension d = 2m + n 
and will appear in C^s with a coefficient of mag nitude g'^'^T^'^. The case n = 1 is an 
exception, because the 1-loop diagram with two external lines is momentum-independent and 
therefore only the operator cfp' is generated at order g"^. Operators of dimension d = 2m + 1 
with the schematic structure V^™0^ have coefficients with magnitude g^T^^"^ for m > 2. The 
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only other exception is 0^, which is generated at tree level and has a coefficient of magnitude 



III. 



Short-distance Coefficients 



The coefficients of the operators in the effective lagrangian (P) must be tuned as functions 
of g, T, and A so that the effective theory reproduces the static correlation functions of the 
full theory at distances 1/T. The parameters can be determined by computing various 
static quantities in the full theory, computing the corresponding quantities in the effective 
theory, and demanding that they match. If the running coupling constant g'^{T) of the full 
theory is small at the scale T, then it is convenient to carry out these calculation using 
perturbation theory in g'^{T). A strict perturbation expansion in g^ is affiicted with infrared 
divergences due to long-range forces mediated by massless particles. These divergences are 
screened at the scale gT, but this screening can only be taken into account by summing 
up infinite sets of higher-order diagrams. This breakdown of perturbation theory does not 
prevent its use as a device for determining the short-distance coefficients in the effective 
lagrangian. As long as we can carry out perturbative calculations in the effective theory 
that make the same incorrect assumptions about the long-distance behavior as perturba- 
tion theory in the full theory, we can match the results and determine the short-distance 
coefficients. 

In the case of the theory, conventional perturbation theory in g"^ corresponds to 
decomposing the lagrangian (^) as £ = Cf^^e + -^int, where 



We will refer to the resulting perturbation theory as a "strict" perturbation expansion in g'^. 
The free part of the lagrangian describes a massless scalar field. A mass will not be generated 
at any finite order in g'^, and the absence of a mass will give rise to infrared divergences that 



■free 




■int 
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become more and more severe as you go to higher and higher orders in g'^. This behavior is 
physically incorrect. It will be clear from the effective theory that a mass m of order gT is 
generated by higher loop diagrams and it provides the screening that cuts off the infrared 
divergences. One way of dealing with the infrared divergences in the full theory is to use 
a reorganization of perturbation theory that incorporates the effects of the mass m in the 



free part of the lagrangian. The simplest possibility |15] is to write the lagrangian (^) as 

C = £free + >Cint, whcre 

Ant = 1^$^ - Im^^'. (9) 

Both g"^ and in the interaction term are treated as perturbation parameters of the same 
order. The mass parameter must be of order g^T"^ in order to avoid large perturbative 
corrections that grow quadratically with T in the high temperature limit. Another possibility 
is to add the mass term to £frcc and subtract it from Ci^t only for the zero-frequency mode 



of $(x, r) Both of these approaches have the drawback that the resulting sums and 

integrals involve two momentum scales T and m, making calculations unnecessarily difficult. 

A simpler approach is to calculate in both the full theory and the effective theory using 
ordinary perturbation theory in g'^, but with an infrared cutoff to regularize the infrared 
divergences. In the full theory, this strict perturbation expansion in g^ is defined by the 
decomposition @. In the effective lagrangian (|^), the coefficients and A are of order g'^ 
and all the coefficients in 6C are of order g'^ or higher. Thus, in the effective theory, the 
strict expansion in g"^ is defined by the decomposition CcS = ('^^efr)free + ('^cfr)int5 where 

(Aff)free = ^(V0)^ 

(£eff)i„t = ^^'(A) 0' + ^A(A) 0^ + SC. (10) 

The expansions in g"^ defined by (P) and ([T0|) both generate infrared divergences that become 
more and more severe in higher orders. But if the parameters in the effective lagrangian are 
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tuned so that the two theories are equivalent at long distances, then the infrared divergences 
in their strict perturbative expansions will also match. Thus, in spite of the fact that the 
strict expansion in gives a physically incorrect treatment of infrared effects, we can use it 
as a device for computing short- distance coefficients. 

Ilia. Coefficient of the unit operator 

In this subsection, we calculate the parameter / in to next-to-next-to- leading order in g"^. 
The parameter / can be interpreted as the coefficient of the unit operator which has been 
omitted from the effective lagrangian (|^). We will determine / by matching calculations of 
log Z in the full theory and in the effective theory. 

We first calculate log^ to next-to-next-to-leading order in g'^ using the perturbation 
expansion for the full theory defined by the decomposition (§) . It is given by the sum of the 
Feynman diagrams in Fig. 1: 



where the sum-integral notation is defined in Appendix A. We regularize both ultraviolet 
and infrared divergences using dimensional regularization in 3 — 2e spatial dimensions, tak- 
ing the momentum scale introduced by dimensional regularization to be A. In (|1^) and 
below, we use the symbol "~" for an equality that holds only in a strict perturbation ex- 
pansion in powers of g'^. Such an equality does not properly take into account the screening 
of infrared divergences at the scale gT , but it can be used for determining short-distance 
coefficients. The sum-integrals appearing in (|11]) are given in Appendix A. To the order that 
is required, renormalization of the coupling constant in the MS scheme is accomplished by 
the substitution 




2 ' 



(11) 



1 + 




(12) 
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After this renormalization, the final result is 



T logZ 
V 
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1^ 
8 167r2 



'A 31 C'(-l) C'(-3) 

31ogy-;^ + — +7 + 4 ^.: .: -2-^^ ^ 



9' 



167r2 



■,(13) 



47rT ■ 15 ■ ' ■ -C(-l) C(-3) 
where 7 is Euler's constant and C,{z) is the Riemann zeta function. The apparent dependence 
of the right side of (|13]) on A is illusory. The renormalization group equation for the coupling 
constant, 

d 9' , ^ 



/i- 
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+ 0{g' 



(14) 



implies that the explicit logarithmic dependence on A in the g'^ term of (|l^) is cancelled by 
the A-dependence of the coupling constant in the term. Thus, up to corrections of order 
g^, we can replace A on the right side of ([T3|) by an arbitrary renormalization scale fi. 
We now consider log Z in the effective theory: 



logZ = -f{A)V + logZ, 



eff- 



(15) 



The partition function for the effective theory is the functional integral in (|^). In the 
diagrammatic expansion for log^efr, the leading terms are given by the diagrams in Fig. 1, 
plus additional diagrams involving mass insertions as in Fig. 2. To match with the strict 
expansion in g'^ for the full theory, we should calculate in the effective theory using the 
perturbation theory defined by the decomposition (p!0D , again using dimensional regulariza- 
tion to regularize both ultraviolet and infrared divergences. This calculation is trivial, since 
massless loop diagrams with no external legs vanish in dimensional regularization due to a 
cancellation between ultraviolet poles in e and infrared poles in e. The result is therefore 

TlogZ 



V 



-fT. 



(16) 



Matching the results (|T3D and ([T6|) , we obtain the coefficient / to order g'^: 
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31 C'(-l) C'(-3) 



AtiT 15 



C(-i) C(-3) 



(17) 



13 



where g'^{fi) is the couphng constant in the MS renormahzation scheme at the scale /z. We 
have used the renormahzation group equation (14) to replace A in (|13|) by an arbitrary scale 
fi associated with the renormahzation of the full 4-dimensional theory. Thus, at this order 
in g"^, the coefficient / does not depend on the ultraviolet cutoff A of the effective theory. 



Illb. Mass parameter 

In this subsection, we calculate the coefficient m^(A) of the 0^/2 term in the effective la- 
grangian (j^) to next-to-leading order in g^. The parameter m(A) can be interpreted as 
the contribution to the screening mass from short distances of order 1/T. The actual 
screening mass rris is defined by the condition that the propagator for spacelike momen- 
tum K = {ko = 0,k) has a pole at = — m^. The physical quantity coincides with 
m^(A) at order (7^, but has corrections of order g^ which arise from the long-distance 
scale l/{gT). In contrast, the mass parameter m^(A) receives contributions only from the 
short- distance scale 1/T, and thus has a perturbative expansion in powers of g'^{T). 

One way to determine m?{h) is to match the propagator for the zero- frequency mode of 
the field $(x, r) in the full theory with the propagator for 0(x) in the effective theory. At 
leading order in g"^, these operators are related as in (|^). This identification is sufficient for 
determining m?{A) to next-to-leading order in g"^. Beyond that order, we must allow for a 
short- distance coefficient multiplying 0(x) in (|^), and we must also allow for the fact that 
</)(x) is only the first term in an operator expansion that contains 0^(x) and other higher 
dimension operators. In general, we must include all operators that are odd under — > —0, 
each multiplied by a short- distance coefficient. In order to determine m^(A) and all the 
necessary short-distance coefficients, we would have to match not only the propagator but 
other 2-point functions as well. 

A simpler way to determine m^(A) is to match the screening mass in the full theory and 
in the effective theory. The screening mass gives the location of the pole in the propagator 
for the zero- frequency mode /q'' (ir$ (x, r) . Denoting the self-energy function for the field 
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$(x, r) at momentum K = (/cq, k) by n(A;o, k), the screening mass is the solution to the 
equation 

P + U{0,k) = at = -ml (18) 

The location of the pole is independent of field redefinitions. Since the operator expansion 
that generalizes (^ can be interpreted as a field redefinition, the screening mass also 
gives the location of the pole in the propagator for the field 0(x). Denoting the self-energy 
for 0(x) by nefj(/c, A), the screening mass must satisfy 

k"^ + m^{A) + U^sik,A) = at k'^ = -m^. (19) 



By matching the expressions for obtained by solving (|T^) and (19), we can determine the 
short-distance parameter m'^{A). 

We will obtain a perturbative expression for the screening mass in the full theory 
by calculating Il{K) to order using the strict perturbation expansion defined by the 
decomposition ([10|) . The self-energy is given by the sum of the Feynman diagrams in Fig. 3: 



u(K) ^ ^ - ^14 ^4 ^ - ^14 1 f20) 

^ 2 TpP^ ATpP^Tp{P^Y Q Tpq P^Q^{P + Q + Kf ^ ' 

We can simplify the equation (|18|) by expanding 11(0, k) as a Taylor expansion around k = 0. 



This can be justified by the fact that the leading order solution to (18) gives a value of k 
that is of order gT, and T^{K) is independent of K at leading order. After setting i^T = in 
the last integral in (pO|), the solution to (18) to order g'^ is trivial: 



2 TpP"^ a Tp P"^ Tp {P'^f Q Tpq P'^Q'^{P + Qf 

It should be emphasized that this perturbative expression does not give a physical value for 
the screening mass, because the sum- integrals are infrared divergent. However, as long as we 
can compute in the effective theory in a way that makes the same incorrect assumptions 
about the long distance behavior, we can match the perturbative expressions to determine 
the short- distance parameter m'^{h). 
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In order to match with the expression (21), we have to calculate the screening mass in 



the effective theory using the strict expansion in defined by the decomposition ([101) . The 
self-energy function T{f,Q{k^ A) in the equation (|19|) for the screening mass has a diagrammatic 
expansion including the diagrams in Fig. 3 and the mass-insertion diagrams in Fig. 4. In 
the full theory, when n(i^) is expanded as a Taylor expansion around if = 0, terms that in 
dimensional regularization scale like fractional powers of k"^ are automatically set to 0. We 
should therefore make the same simplifications in the effective theory. But then all the loop 
diagrams in Fig. 3 and Fig. 4 vanish, since the external momentum k provides the only mass 
scale in the integrals. The self-energy function reduces to 

neff(A;,A) ^ 5m\ (22) 

where 5m? is the mass counterterm that contains the poles in e that are associated with 
mass renormalization. The solution to the equation (|19|) for the screening mass is therefore 
trivial: 

ml ^ m^(A) + 6m^. (23) 

From (P5|), we see that the screening mass in this unphysical perturbation expansion is just 
the bare mass. Comparing (PT|) and (PSf), we find that m^(A) is given by 



2 fpP^ AfpP'^fp{P^Y Q fpQ P^Q^{P + Qy 
where the sum-integrals are to be evaluated using dimensional regularization of both ultra- 
violet and infrared divergences. The sum-integrals in (^) are given in Appendix A. After 
renormalization of the coupling constant g, there remains a pole in e which must be cancelled 
by the mass counterterm Sm"^. The mass counterterm is thereby determined to be 

Sm^ = —i (25) 

24(167r2) e ^ ^ 

Our final expression for the short- distance mass parameter m^(A) is 
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m\A) = ^Af^)T' 1 + 



-31og^ + 41og^ + 2-, + 2-^^_^^ 
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where g^{fi) is the MS couphng constant. We have used the renormahzation group equation 



(H) to change the renormahzation scale of the full theory from A to /x. The remaining 
logarithm of A in (26) reveals that m^(A) depends explicitly on the factorization scale A at 
order g'^. This A-dependence is necessary to cancel logarithmic ultraviolet divergences from 
loop integrals in the effective theory. 



IIIc. Coupling constants 

For the calculations in this paper, we require the coupling constant A of the 0^ interaction in 
the effective theory only to leading order in g^. At this order, we can simply read A off from 
the lagrangian of the full theory. Substituting $(r, x) — > \/T(j){yi) in and comparing 
Jq drC with Ccs in (H), we find that, to leading order in g"^, 

A = g'T. (27) 

There is no dependence on the factorization scale A at this order. The coupling constant A 
could be calculated to higher order in g"^ by matching 4-point correlation functions in the 
full theory and in the effective theory. Beyond next-to-leading order in g"^, the matching is 
complicated by the breakdown of the simple relation (^) between 0(x) and the fundamental 
field. A more convenient quantity for matching beyond leading order is the on-shell scattering 
amplitude defined by the residue of the 4-point function at the poles of the propagators of 
the four external lines. Like the screening mass, this scattering amplitude is invariant under 
field redefinitions. 

The only other coefficient in the effective lagrangian that is known is the coefficient of 
0^. It has been computed by Landsman [T^, and its value is 15C(3)(7®/(1287r'^). It first 
contributes to the free energy density at order g^ and to the square of the screening mass at 
order g^. 
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IV. Calculations in the Effective Theory 

In this section, we calculate physical quantities in the effective theory using perturbation 
theory. We calculate the free energy to order , reproducing a recent result by Parwani and 
Singh 0, and we obtain a new result for the screening mass to order g^. 

IVa. Free Energy to Order 

The free energy density F{T) is defined by 

Z{T) = e-^^(^)^. (28) 
Comparing with the equivalent expression for the partition function (0), we obtain 

F{T) = T /(A) - T , (29) 

where the partition function for the effective theory is 

^cff = Vcj) exp 1^- J d^x Cft^ . (30) 



The strict perturbation expansion for log Z^s corresponding to the decomposition (|ToD of Ces 
contains infrared divergences. These divergences were not a problem in the matching calcu- 
lations of Section III, since identical infrared divergences appeared in the strict perturbation 
expansion for the full theory. However, if we wish to actually calculate the free energy, we 
must incorporate the physical effects that cut off the infrared divergences into the free part 
of the lagrangian. The necessary infrared cutoff is provided by the 0^ term in the effective 
lagrangian. We therefore make the following decomposition of Ces into free and interacting 
parts: 

M^ee = ^(V0)' + ^m2(A)0^ 

i^eS^nt = ^A(A) 0^ + 6C. (31) 

From the matching calculations in section III, is of order g'^T'^ and A is of order g^T. 
Since the only momentum scale in (£efr)free is m, any powers of T in the coefficient of an 

18 



operator will be compensated by powers of m. Thus the effective expansion parameter for 
the 0^ perturbation is A/m, which is of order g. The next most important perturbation after 
0^ is the dimension-4 operator (0V0)^, for which the effective expansion parameter is of 
order g'^m/T or, equivalently, of order g^. Similarly, the effective expansion parameter for 
the dimension-3 operator 0^ is g^. Thus to calculate logZcs to next-to-next-to-leading order 
in g, we need only consider the 0^ perturbations. 

The contributions to logZ^s of orders g^, g'^, and g^ are given by the sum of the 1-loop, 
2-loop, and 3-loop diagrams in Fig. 1 and the first diagram in Fig. 2: 



logZeff 1 /", / 2 , 2^ ^ 



2 



AV /• 1 V r 1 

+ T77 



16 \Jp + m? ) Jp {p^ + m?Y 

+ I . (32) 

48 Jpqr {p^ + m?) (g^ + m?){r^ + m?) [(p + q + r)^ + m?] 

The only effect of the diagram with the mass counterterm in Fig. 2 is to replace m? in the 
1-loop diagram by mo = m? + 5m?, where the mass counterterm is 



With the identification A = g'^T , this expression is identical to the expression in (pSf). When 
/log(p^ + "^o) is expanded in powers of A using (^), the 1/e term cancels against a pole 
from the last integral in (0), but it also gives rise to finite contributions from the expansion 
of the integral /log(p^ -|- m^) to order e. The integrals in (|32D are given in Appendix B. 
Adding up the diagrams, we obtain 



logZeff 1 3... 1 A 2 1 
— -m (A) — m — 
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2 



m . (34) 



V 127r ' ' 87rl67r 127r 

Our final result for the free energy ( PP| ) to order g^ is the sum of two terms that represent 
the contributions from the momentum scales T and gT, respectively: 

F(r) = Fi(r) + F2(T). (35) 
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The first term Fi(T) = fT is the contribution to the free energy from single-particle effects 
involving the short-distance scale 1/T only. This contribution can be expressed as a power 
series in g'^(2TTT), with the first three terms given by ([T7|): 



Fi(T) 



lg\27rT) 



_ _ - : 

9 1 10 8 167r2 



+ 



167r2 



(36) 



The choice fi = 2ttT for the renormalization scale of the effective theory will be justified 
in Section V. The second term F2{T) = —TlogZcs/V in the free energy (^) takes into 
account collective effects of the particles involving the long-distance scale l/{gT). It can be 
expressed as a perturbation series in A and in the other coupling constants of the effective 
theory. The first three terms are given by (P^: 

2 \ 



F2{T) 



^-m3(A)r|l ^ ^ 



A 9 
41og— + - - 41og2 
2m 2 
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IGnm 



(37) 
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where A = g'^T and the mass parameter m(A) is given by ( p6l) with fi = 2ttT. Note that, 
after substituting the expression (|26|) for the mass parameter m^(A), the A-dependence in 
(^) cancels to next-to-next-to-leading order in g. Adding the short-distance contribution in 
(^) and expanding in powers of g{27cT), the free energy reduces to 



F{T) 



IT 



9 \ 10 8 V 47r ; VQ\ 4:7V J 



1 

+ 7^ 



59 , C'(-l) C'(-3) 
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(38) 



The coefficient of g'^ was first calculated by Frenkel, Saa, and Taylor |]l[, up to an error 
that was corrected by Arnold and Zhai 0. The order g^ term was recently calculated by 
Parwani and Singh . Our result agrees with theirs after taking into account the difference 
in the definition of the coupling constant. There is a loss of accuracy when we make a strict 
expansion in powers of g as in (|38|) . In Section V, we will give a more accurate expression 
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for the free energy which sums up leading logarithms of g from higher orders of perturbation 
theory. 



IVb. Screening Mass to Order 

The screening mass rUg describes the long-distance behavior of the potential produced by 
the exchange of a particle with spacelike momentum. The potential falls exponentially like 
g-rrisH large R. The screening mass is a long distance quantity, so it should be calculable 
using the effective field theory. In this section, we use the effective field theory to calculate 
rris to order g"^. 

The screening mass m^, which gives the location of the pole in the propagator of the 
effective theory, is the solution to the equation (|T^. At leading order in A/m, the solution 
is simply = m(A). The self energy function Iles{k,A) is given to next-to- leading order 
in A by the Feynman diagrams in Figs. 3 and 4: 

TT A^ A /• 1 X' r 1 r 1 



2 Jp p'^ + m? 4 Jp p^ + m? Jp {p^ + m?Y 

6 Jpq {p^ + m?) (g^ + w?) [(p + q + k)^ + rn?] ~^ ^ ^ 

The mass counterterm 5m?, which is given in (^), cancels an ultraviolet pole in e in the 
integral over p and q. This is the only integral in ( PU] ) that depends on k. The self-consistent 
solution to (JT^) to next-to- leading order in A/m is obtained by evaluating the integral at the 
point k = im. The resulting expression for the screening mass is 

X^ r 1 r 1 



ml 



2 Jpp^ + ' 



m^ 4 Jp p'^ + m^ Jp {p^ + m^y 
X' r 1 



+ 6m\ (40) 

k=im 



6 Jpq {p^ + m'^){q^ + m^) [(p -|- q -|- k)^ -|- 
Note that the screening mass is not identical to the value of the inverse propagator at 
momentum, which would be given by (^0[) with the last integral evaluated at = 0. Unlike 
the screening mass, the mass defined by the inverse propagator at A; = is not invariant 
under field redefinitions. 
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The integrals in ( ^0]) are given in Appendix B. To next-to-next-to- leading order in A/m, 
the screening mass is 

2 \ 
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IGvrm 



41og^ + 3-81og2 
2m 
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(41) 



Note that, when we substitute p6| ) for m^(A), the A-dependence cancels to next-to-next- 
to- leading order in g. This verifies that the screening mass is independent of the arbitrary 



factorization scale A to this order. Setting /i = 27rT in (j26|), substituting it into (^T]), and 
expanding in powers of g{27!-T), we obtain 



ml = i-/(27rT)T2 Jl- 
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(42) 



The term of order in (^5|) was first calculated by Dolan and Jackiw [rB|. The correction 
of order g^ in the expression for the screening mass is a new result. The screening mass 
differs at order g^ from the quasiparticle mass, which was calculated to order g'^ by Parwani 
rS|. The quasiparticle mass is defined by the pole in the energy for the propagator at zero 



3-momentum. It gives the energy of the single-particle excitations in the plasma, while the 
screening mass gives the range of the force mediated by the exchange of a particle in the 
plasma. There is a loss of accuracy in making a strict expansion of the screening mass in 
powers of g as in (^). In Section V, we will give a more accurate expression for the screening 
mass which resums leading logarithms of g from higher orders of perturbation theory. 



V. Summation of Leading Logarithms 

The coefficients in the effective lagrangian (^ depend on not only on the coupling constant 
g'^ and the temperature T, but also on two arbitrary momentum scales: the renormalization 
scale /i of the full theory and the ultraviolet cutoff A of the effective theory. In this section, 
we exploit these arbitrary scales to sum up the leading logarithms from higher orders in 
perturbation theory. We first discuss the choice of the renormalization scale /x. We then 
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present the evolution equations that describe the dependence of the short- distance coefficients 
on the scale A. We then show how the solutions of the evolution equations can be used to 
sum up leading logarithms of T / (gT) in physical quantities. 

Va. Renormalization scale 

The short-distance coefficients in the effective lagrangian would be independent of the ar- 
bitrary renormalization scale /i of the original theory if they were calculated to all orders 
in g^. A dependence on /i appears however when a coefficient is calculated only to a finite 
order in g^. At leading nontrivial order in g"^, the scale /i appears to be completely arbitrary. 
The resulting ambiguity can be decreased by a next-to-leading order calculation, because 
an inappropriate choice of fi will result in unnecessarily large perturbative corrections. The 
short- distance coefficients / and m^(A) were both calculated to order g'^ in Section III, and 
we can use those results to discuss the appropriate scale for /i in the g'^ terms. 

The transcendental constants appearing in the expressions (p!7[) and (|26|) have the nu- 
merical values 

C'(—l) (T'f— 3) 

7 = 0.57722, — f = 1.98505, ^ — f = 0.64543. (43) 

^ C(-i) C(-3) ^ ^ 



The short-distance coefficients then reduce to 



/ = --TU1 



90 4 167r2 



1 + ( -7.2 - 3 log ^ ^ ^ 
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-^(A) = Y/if^)T' 



2'kTJ 167r2_ 

1+ (4.7-31og^ + 41og^)^ 



(44) 
(45) 



where the coupling constant is the MS constant at the scale fi. This couphng constant is 
defined so that it is the appropriate coupling constant for particles whose Euclidean invariant 
mass-squared k^ + k"^ is approximately /x^. The parameters (44) and ( ^5]) are coefficients in an 
effective lagrangian obtained by integrating out modes with nonzero Matsubara frequencies 
kf) = 2n'n'T, n = 1,2,3,.... Since these modes have invariant masses that are equal to or 
greater than 27rT, we expect n = 2t[T to be an appropriate choice for the renormalization 
scale. Choosing = A = 27iT, the coefficients of 5'^/(167r^) in square brackets in (|^) and 
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in (|5) are —7.2 and 4.7, respectively, which are reasonably small. The more naive choice 
H = A = T gives even smaller coefficients —1.7 and 2.9. We prefer the physically motivated 
choice, and we therefore set the renormalization scale to /i = 27rT throughout the remainder 
of this paper. 



Vb. Factorization scale 

The matching calculations described in Section III give the short-distance coefficients in the 
effective lagrangian C^s as perturbation series in g^{2TTT) with coefficients that are polynomi- 
als in log(27rT/A). To avoid unnecessarily large coefficients in these perturbative expansions, 
we must choose A of order 27rT. Once the short- distance coefficients are known, physical 
quantities can be calculated in the effective theory using a perturbation expansion in A/m 
and in other dimensionless parameters obtained by dividing the short- distance parameters 
by appropriate powers of m. The coefficients in the resulting perturbation expansions for 
physical quantities contain logarithms of A/m. To avoid unnecessarily large coefficients in 
these perturbation expansions, it is necessary to carry out the perturbative calculations in 
the effective theory using short-distance parameters that are evaluated at a scale A of order 
m. The parameters calculated at the original scale A = 27rT must therefore be evolved 
down to the scale A = m before they can be used in these perturbative calculations. The 
A-dependence of these parameters is described by "renormalization group equations" or 
"evolution equations" . 

The effective lagrangian (|^) can be expressed as a sum over over all local operators that 
respect the symmetries of the theory: 

/(A) + C,s = ECn(A)a, (46) 

n 

where we have included the unit operator as one of the operators On- The coefficients C„ 
are the generalized coupling constants of the effective theory. Because of ultraviolet diver- 
gences, the effective theory must be regularized with an ultraviolet cutoff A. The ultraviolet 
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divergences in the effective theory include power ultraviolet divergences proportional to A^, 
p = 1,2, . . ., and logarithmic divergences proportional to log(A/m). As discussed in Section 
II, the power divergences are artifacts of the regularization scheme and have no physical con- 
tent. If they are not removed as part of the regularization procedure, they must be cancelled 
by power divergences in the coupling constants C„. In contrast, the logarithmic ultraviolet 
divergences are directly related to logarithms of T in the full theory, and therefore repre- 
sent real physical effects. This difference justifies treating power ultraviolet divergences and 
the logarithmic ultraviolet divergences differently. It is convenient to use a regularization 
procedure for the effective theory in which power ultraviolet divergences are subtracted and 
logarithmically ultraviolet divergent integrals are cut off at the scale A. These logarithmic 
divergences are then the only ones that must be cancelled by the A-dependence of the cou- 
pling constants C„. The dimensions of a coupling constant can then only be taken up by 
powers of the temperature T. The coupling constant C„ must be proportional to T^~'^^ 
where dn is the scaling dimension of the corresponding operator On- The dimensionless fac- 
tor multipling T^-^^ ^j^g coupling constant C„ can be computed as a perturbation series in 
g'^{T), with coefficients that are polynomials in log(T/A). The dependence on A is governed 
by a "renormalization group equation" or "evolution equation" of the form 

A^a(A) = PniC), (47) 

where the beta function /3„ has a power series expansion in the coupling constants C^. 
These equations follow from the condition that physical quantities must be independent of 
the arbitrary scale A. 

Since C„ is proportional to T^"*^" , every term in the expansion of its beta function must 
be proportional to T^~'^^. In particular, a term like Cm^Cm2 ■ ■ ■ Cm^. can appear only if the 
dimensions d^. of the corresponding operators O^. satisfy 

k 

=3-4. (48) 

i=l 

This condition is very restrictive, particularly if the effective field theory is truncated to those 
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terms that are given explicitly in (^. It implies that the beta function for the coefficient / 
of the unit operator can only have two terms proportional to A'^ and m^A. The beta function 
for must be proportional to A^ and the beta function for A must vanish to all orders in A. 
These restrictions reflect the super-renormalizability of this truncated effective theory, which 
implies that there are only a finite number of independent ultraviolet-divergent subdiagrams. 

The evolution equations for / and are calculated in Appendix C. They follow from 
the condition that the free energy and the screening mass must be independent of A. The 
evolution equations are 

A±X ^ 0. (51) 

The allowed term m^A does not appear in the beta function for /. Using A = g^T, we see 
that the beta function for in (|50|) is consistent with the explicit calculation to order 
in (1^). The evolution equations (^91), (0), and (|5l|) are correct to all orders in A, but 
they receive corrections involving the coefficients of higher dimension operators like (p^. The 
coefficient of (p^ in the effective lagrangian is of order g^. It gives corrections to the right 
sides of the evolution equations (^), (^Of), and (|5lD that are of order g^^X^, g^X"^, and g^X, 
respectively. 



To the accuracy given in (^91), (^OD, and (|5TD , the solutions to the evolution equations 
are trivial. The coupling constant A does not evolve with A, and the solutions for / and 
are simply linear in log(A): 

/(A) = /(27rT) - ^ (^A^ iog(A/27rT), (52) 

m^A) = m'{27rT) + ^ (^A^ log(A/27rT). (53) 

Note that, since the perturbative expression (^) is linear in log(A), it already satisfies the 
evolution equation (|50D. 
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Vc. Resumming Logarithms of g 

The evolution equations for the short- distance coefficients can be used to sum up leading 
logarithms of T/{gT) in physical quantities, such as the free energy and the screening mass. 
We must first choose a value for the scale A which will avoid unnecessarily large coefficients in 
the perturbation expansions of the effective theory. A reasonable choice is the screening mass 
m<j, since this is the minimum invariant mass for a particle in the 3-dimensional Euclidean 
theory. At leading order in g, the screening mass is simply = gT/\/2A^ so we choose 
A = gT/^/2A. 

An expansion for the screening mass in powers of g is given in (^). This expression is 
only accurate up to corrections of order g^ log (7. A more accurate expression for the screening 
mass can be obtained by not expanding out the mass parameter in (|4l]) : 

where g = g{2TTT) is the coupling constant at the scale 2itT and is the mass parameter 
in dH) with /i = 27tT and A = gT/V2A: 



41og^^ + 2-log2-7 + 2-^*' ^' 
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The expression (^4D is accurate up to corrections of order g^ . It also sums up all the "leading 
logarithms" of order g'^^^^ log" (7, n = 2, 3, . . .. These terms are generated by expanding out 
the term proportional to g'^Tm in (|5^) using (p5|). 

An expansion for the free energy in powers of g is given in (^Sj). This expression is 
accurate up to corrections of order g^ log (7. A more accurate result can be obtained by using 
the solution to the renormalization group equation (^91) for / in the short-distance part and 
by not expanding out the mass parameter m in the long-distance part 
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where g = g{2'KT), and m is given by (0). The first two terms on the right side of (|56|) are the 
short distance contribution /(A)T, with the factorization scale evaluated at A = gT/\/2A. 
Its expansion to order g'^ is given by Fi{T) in (|36D, and the term of order g^logg comes 
from the solution (0) to the renormalization group equation for /. There is also another 
contribution of order g^logg that comes from expanding out the g'^im?T'^ term in (|56|) using 
(^). Having included both of these g^logg terms, the expression (56) for the free energy is 
accurate up to corrections of order g^ . It also sums up all the "leading logarithms" of the 
form ^f^""^^ log" ^f, n = 2, 3, . . ., which are obtained by expanding out the nv'T term in (|56|). 



VI. Conclusions 

We have developed an effective-field-theory approach for calculating the thermodynamic 
properties of a field theory in the high temperature limit. The effective field theory is the 3- 
dimensional theory obtained by dimensional reduction to the bosonic zero-frequency modes. 
The short-distance coefficients in the effective lagrangian are computed by straightforward 
perturbative calculations in the full theory without any resummation. Thermodynamic 
quantities are then calculated using a perturbation expansion in the effective theory which 
incorporates the effects of screening. In each of these two steps, the calculations involve 
only a single mass scale, which greatly simplifies the sums and integrals that need to be 
evaluated. The short- distance coefficients satisfy renormalization group equations which can 
be used to improve the perturbation expansion by summing up leading logarithms of T / {gT). 
The power of this method was demonstrated by carrying out two calculations in massless 

theory beyond the highest orders that were previously available. The free energy was 
calculated to order g^logg and the screening mass was calculated to order g^logg. 

The method that we have used to calculate the free energy of a scalar field theory can also 
be applied to nonabelian gauge theories, such as QCD. The free energy for QCD has been 
calculated with an error of order g^logg using more conventional resummation methods [Q. 
Using our effective field theory approach, it should be straightforward to decrease the error 
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to order g^. The effective theory obtained by integrating out the scale T is a 3- dimensional 
gauge theory with a scalar field in the adjoint representation. There are two short distance 
coefficients that must be calculated beyond leading order in g"^ in order to compute the free 
energy to order g^. The coefficient / of the unit operator is required to next-to- next-to- 
leading order in g"^, but this is already known [Q. The electric mass parameter m^j must be 
calculated to next-to- leading order in g"^. The error can be decreased further by calculating 
the renormalization group equations for / amd m^i. Their solutions can be used to sum up 
terms of the form g^log^ g and g^logg, thereby reducing the error to order g^. This is the 
maximal accuracy that can be achieved using purely diagrammatic methods. At order g^, 
there is a contribution to the free energy from the momentum scale g^T that can only be 
calculated using lattice simulations of 3-dimensional QCD ||^. 

A similar effective field theory approach was recently developed by Farakos, Kajantie, 
Rummukainen, and Shaposhnikov to study the electroweak phase transition [|T7|. They 
integrated out the scale T to get a 3-dimensional effective theory, and exploited the renor- 
malization group equations of the effective theory to take into account logarithms of T/ (gT), 
just as we have done in this paper. They also integrating out the scale gT to obtain a second 
effective field theory that must be treated numerically. This second step is necessary in gauge 
theories since there is no screening of magnetostatic forces at the scale gT. This strategy 
has also been used to determine the asymptotic behavior of the correlator of Polyakov loops 
1^ and to solve the problem of calculating the magnetostatic contribution to the free energy 
of a nonabelian gauge theory [Q. 

An effective field theory approach has also been applied recently to the massless $^ 
theory by Marini and Burgess |T3. They used a momentum cutoff as a regulator and their 
short-distance coefficients therefore contain power ultraviolet divergences that serve simply 
to cancel the power ultraviolet divergences from loop integrals. These power divergences 
greatly complicate the renormalization group equations for the short- distance coefficients. 
As we have emphasized, the power divergences are artifacts of the regulator and might as 



29 



well be subtracted as part of the regularization scheme. This not only greatly streamlines 
explicit calculations, but it also makes the conceptual framework more transparent. 

Effective field theories obtained by dimensional reduction have provided great insight 
into the qualitative behavior of field theories in the high temperature limit. In this paper, 
we have shown how they can also be used as a practical tool for explicit calculations. By 
using effective field theory to separate the important momentum scales T and gT (and g^T 
if necessary) , perturbative calculations can be organized into steps that involve only a single 
momentum scale at a time. By exploiting the renormalization group structure of the effective 
theory to sum up leading logarithms of g, potentially large coefficients in the perturbative 
expansion can be brought under control. We have exhibited the power of our effective-field- 
theory method by carrying out pioneering calculations in massless $^ theory. This method 
has many other exciting applications, especially in unravelling the complexities of nonabelian 
gauge theories at high temperature. 
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A Sum- integrals in the Full Theory 

In the imaginary-time formalism for thermal field theory, a boson has Euclidean 4- 
momentum P — {po,p), with — Po + P^- The Euclidean energy po has discrete values: 
Po — 2TrnT, where n is an integer. Loop diagrams involve sums over po ^-nd integrals over p. 
It is convenient to introduce a concise notation for these sums and integrals. If dimensional 
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regularization is used to regularize ultraviolet or infrared divergences, the definition is 
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where 3 — 2e is the dimension of space and fi is an arbitrary momentum scale. The factor 
[e'^ /Any is introduced so that, after minimal subtraction of the poles in e due to ultraviolet 
divergences, fi coincides with the renormalization scale of the MS renormalization scheme. 

The sum-integrals required to calculate the coefficient /(A) to next-to- leading order in 
g"^ can be found in Ref. We reproduce them here for convenience: 
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where 7 is Euler's constant and ({z) is Riemann's zeta function. 
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B Integrals in the Effective Theory 

The effective theory for the scale g'^T is an Euclidean field theory in 3 space dimensions. 
Loop diagrams involve integrals over 3- momenta. It's convenient to introduce the notation 
/p for these integrals. If dimensional regularization in 3 — 2e dimensions is used to regularize 
ultraviolet divergences, we use the integration measure 



p y 4n J J (27r)3-2^ ■ ^^'^^ 
If renormalization is accomplished by the minimal substraction of poles in e, then /i is the 
renormalization scale in the MS scheme. 
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The integrals that are required to calculate the free energy to order and the screening 
mass to order g'^ are 

,3 



/ log(j9^ + 

■J V 



pP^ + mi? 
1 



m 

4:71 
1 



Jp {p^ + m?Y Svrm 
1 1 1 

pq + IV? (f + vn? (p + q)2 + vn? 



2rr 

l + .(21og£)+0(e=) 



67r 



- + 41og-^ + 2 + 41og2-41og3 + 0(e) 
e 2m 



(B.2) 
(B.3) 
(B.4) 

(B.5) 
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fe=im 

- + 4 log — + 6 - 8 log 2 + 0(e) 
e 2m 
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1 



1 



pqr + m^ + m? + (p + q + r)^ + 

m 



(47r) 



- + 61og — + 8 - 41og2 + 0(e) 
e 2m 
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(B.7) 



The integral (|B.2|) is standard. The remaining integrals are most easily evaluated by 
going to coordinate space. The Fourier transform of the propagator l/(fc^ + m^) in 3 — 2e 
dimensions defines a potential V{R): 

1 



V{R) ^ / e*-^ -— ^ 
Jk + rw^ 



(B.8) 



It can be expressed in terms of the modified Bessel function Ky{z): 

1 /m\V2-€ 



V{R) 



^ 47r ; (27r)3/2- 
At e = 0, this reduces to the familiar Coulomb potential 



(B.9) 



Vo{R) 



-mR 



AttR 



(B.IO) 



For small R, the potential V{R) can be expressed as the sum of two Laurent expansions in 
R^, a singular one beginning with an term and a regular one begining with an 
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term: 



V{R) 



e) 1 



r(l) 4vr 



R 



-l+2e 



— (e^/x ) „ . 1 . —m 



r(-l) 47r- 



1 + 
1 + 



2(1 + 2e) 
m^R^ 



+ 0{m^R^) 
+ Oim^R^' 



(B.ii: 



2(3 - 2e) 

The integrals ( p.3|) and ( |B.4|) are related to the potential at i? = 0. By the rules of 
dimensional regularization, they can be evaluated by taking e large enough that the i?"^"*"^"^ 
term vanishes as i? — 0, and then analytically continuing back to small values of e. The two 
integrals are 

1 



p p2 + m? 



r(o) 



/p (p2 + m^)^ 
The integrals 



2m (im 



(eV)^ 



r(l + e) 1 



i?=0 



r(^) 



Stt 



-m 



-l-2e 



:b.i2) 
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and ( [B.Tp require a little more effort. The integral ( [B.7| ) can be written 
111 1 



pgr + m? + m? + (p + q + i")^ + 



y^(/2) 



where is defined by 



R 



47r 



(B.14) 



:B.15) 



From the i? — region of the integral (p3.14| ), there is a linear ultraviolet divergence, which 
is removed by dimensional regularization, and a logarithmic divergence, which appears as 
a pole in e. We evaluate the integral by splitting the radial integration into two regions, 
< -R < r and r < R < oo. Since the ultraviolet divergences come only from the region 
i? — > 0, we can set e = in the region r < R < oo. Thus the integral can be written 



V\R) 



- An dRR'-''V\R)+ An dRR'V^{R) . (B.16) 
e) JO 



iR ■ V 4 ; r(^ , 

By choosing r ^ 1/m, we can evaluate the first integral on the right side of ( [B.16| ) by using 
the small- i? expansion ( |B.11| ) for V{R). Dropping all terms that vanish as r — 0, we get 



r 



ri^An / dRR^-^^V\R) 
^-e) Jo 



(47r)3 



1 /I rV ^ 

m - + 2 log — ^ + 47 + 4 

r V e 2m , 



+ 0(e). (B.17) 
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The integral over the region r < i? < oo is easily evaluated using integration by parts. 
Dropping all terms that vanish as r — >^ 0, we get 



dRR^ V\R) 



{Any 



— h 4m (log 4mr + 7 — 1) 



:b.i81 



Note that the 1/r and logr terms cancel between ( |B.17|) and ( B.18 ). Inserting ( |B.17| ) 
and (IHAHD into (IKTeD , we obtain 

The integral ( [B.6|) can be evaluated in a similar way to ( |B.7| ). It can be written 

11 1 



(B.19) 



ipq p2 + g2 _|_ ,7^2 (p + q + k)^ + Jr 

Again we split the radial integration into two regions and set e = in the region r < R < 00. 
After evaluating the angular integrals, we obtain 

^ (277)3/2 
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^, , ^ r dRR^/^~'Jy2-e{kR)V^R) 
2k J y/k Jo 

+ — / dRRsm{kR) V^{R) , 
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where Ju{z) is an ordinary Bessel function. We evaluate the first integral using the small- i? 
expansion ( p.ll|) for V{R) and the small- i? expansion for the Bessel function 



Jl/2-e{kR) 



Til-e) \2 



-kR 



1/2-e 



l + 0(P/^2^ 



:b.2i] 



Dropping terms that vanish as r ^ 0, the first integral in ( p.20|) is 



2k 



^ r dR R^'^-' Ji/2-e{kR) V'\R) 
\jk h 

^ ^- + 41og/xr + 2 + 47') + 0(e). 
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This integral is independent of k. In the second integral in ( |B.2(j| ), we have to set k = im: 



k J, 



dRRsm{kR) V^{R) 



k=im 



[— log 2mr + 1 — 2 log 2 — 7) 



(B.23) 



Adding ([B.22|) and (pj.23|), the logarithms of r cancel and we obtain 
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C Evolution Equations for Coefficients 

In this appendix, we calculate the evolution equations for the short-distance coefficients 
/ and in the effective lagrangian. These equations follow from the condition that phys- 
ical quantities must be independent of the arbitrary renormalization scale A of the effec- 
tive theory. The A-dependence of the parameters in the effective lagrangian must cancel 
the A-dependence from loop integrals. If power ultraviolet divergences are subtracted as 
part of the regularization scheme, then the A-dependence comes only from logarithmically- 
divergent loop integrals. One-loop diagrams in a 3-dimensional field theory never give loga- 
rithmic ultraviolet divergences. After averaging over angles, such an integral has the behavior 
/ at large p. This gives a power divergence for n = 1 or less and is convergent for 

n = 2 or greater. Thus logarithmic divergences only arise in diagrams with 2 or more loops. 
If the number of loops is odd, logarithmic divergences only arise from subdiagrams with an 
even number of loops. Thus the evolution equations are completely determined by diagrams 
with an even number of loops. 

The renormalization group equations for can be determined from the condition that 
the screening mass rris is independent of A. Since only the ultraviolet region of loop inte- 
grals is relevant for determining the evolution equations, we can use the expression for the 
screening mass that is obtained from the strict perturbation expansion in g"^ defined by the 
decomposition (0). The expression for the screening mass to order g"^ is 

2 _ 2 A /• 1 _ /• 1 /• 1 _ A^ /• 1 _ r 1 . 

m, fti m + 2 ^ " T ip ^ i ~ y Jpg pV(p + q)' 2~ ip W ^ ' 

(C.l) 

The only logarithmic divergence comes from the 2-loop integral over p and q, which comes 
from the last diagram in Fig. 3. Thus the condition that is independent of A reduces to 



The ultraviolet divergence in the integral on the right side of (|C2|) is the same as in 
The derivative is with respect to the scale A associated with the ultraviolet cutoff, and 

35 



must be taken with the infrared cutoff fixed. The integral in ( U.2 ) vanishes in dimensional 



regularization only if we use the same regularization scale /i for ultraviolet and infrared 
divergences. If we use a different scale A for the regularization of ultraviolet divergences, the 
integral is 

41og-l . (C.3) 



pgpV(p + q)2 (87r)2 



The subscripts UV and IR indicate whether the pole in e is of ultraviolet or infrared origin. 
Inserting ( |C.3| ) into (|C.2D , we obtain the evolution equation 

A—m^ = - . (C.4) 

This result is accurate to all orders in A and to leading order in the coefficients of higher 
dimension operators. 

The evolution equation for / can be determined from the condition that the free energy, 
or equivalently, the logarithm of the partition function given in (15), is independent of A. We 



need only consider diagrams for log Z^s which have an even number of loops, and they can be 



calculated using the strict perturbation expansion in g"^ defined by the decomposition ([To|) . 
The two-loop diagram for log^es in Fig- 1 has no logarithmic ultraviolet divergence. The 
four-loop diagrams that have logarithmic ultraviolet divergences are shown in Fig. 5. The first 
diagram in Fig. 5 has two- loop sub diagrams that are logarithmically divergent. The resulting 
A-dependence is cancelled by the A-dependence of m? in the coefficient of the second diagram 
of Fig. 2. The second diagram in Fig. 5 has no logarithmically divergent subdiagrams, but 
it has an overall logarithmic divergence. The A-dependence of this contribution to the free 
energy can only be cancelled by that of the short-distance coefficient /: 
l^d d r 1 

dA^ ~ "48 dA 4,,,,3 g?(p + qi)2gi(p + q2)^g|(p + qs)' ' ^ ■ ^ 

After combining pairs of propagators using the Feynman parameter trick, the integrals over 
qi, q2, and can be evaluated analytically using ( p.l3|) . The result is 

1 



m<?2<73 Qiip + qO^gKp + q2)^g3(p + qs)^ 
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(C.6) 



The integral over p vanishes in dimensional regularization if the same scale ^ is used in 
the regularization of ultraviolet and infrared divergences. If a different scale A is used for 
ultraviolet divergences, the value of the integral is 



2^2/ 



32(167r) 



-) -(-) +81og- 



Inserting this result into ( |C.6| ), we obtain the evolution equation 

dk^ 12 \ 167r/ 



(C.7) 



(Ci 



This result is accurate to all orders in A and to leading order in the coefficients of higher 
dimension operators. 
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Figure Captions 

1 . Feynman diagrams for the logarithm of the partition function in the full theory and in 
the effective theory. 

2. Additional Feynman diagrams for the logarithm of the partition function in the effective 
theory. 

3. Feynman diagrams for the self-energy in the full theory and in the effective theory. 

4. Additional Feynman diagram for the self-energy in the effective theory. 

5. Four- loop diagrams for the logarithm of the partition function of the effective theory 
which have logarithmic ultraviolet divergences. 
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Figure 1 




Figure 2 
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Figure 4 
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